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A COMPUTER PROGRAM FOR THE THERMAL ANALYSIS OF THE 
FIRE ENDURANCE OF CONSTRUCTION WALLS 


Prancis.c.JW. Fung 


Abstract 


A general one-dimensional transient heat and mass transfer 
numerical program has been developed for composite building con- 
structions. Since typical building constructions consist of a 
series of composite layers and intermediate air layers, transient 
heat transfer is modeled by conduction through solids and by 
radiation and convection through air spaces. In addition the 
program has built-in features for ease of application to 
building constructions where various combinations of solid-to- 
solid and solid-to-air interfaces are encountered. The com- 
plete Fortran language program as used on the NBS Univac 1108 
computer is given. A discussion of the program and instruction 
for its use are facilitated by the aid of examples. Numerical 
solutions using the present program compare favorably with experi- 


mental data in standard fire endurance tests. 


Key words: Composite building constructions; convection; 
experimental data; fire endurance; heat generation or absorption; 
numerical solutions; radiation; one-dimensional; thermal con- 


duction; transient heat transfer. 


1. INTRODUCTION 


The fire performance of building constructions is generally evaluated 
by a large scale laboratory fire endurance test (ASTM E 119) in which one 
surface is exposed to furnace gases with a prescribed increasing temperature 
history simulating the burnout of combustibles. The fire endurance rating 
of the construction is the time period during which it withstands the fire 
exposure without (a) structural failure, (b) the development of cracks 
through which flames or hot gases can pass, or (c) the temperature rise on the 
unexposed surface exceeding, a prescribed limit (250 °F average rise or 325 °F 
rise at a single point). Where the failure criterion is due to heat trans- 


mission alone, heat transfer analysis can provide a means for prediction and 


design. 


Since typical building constructions consist of a series of composite 
layers and intermediate air layers, transient heat transfer is modeled by 
conduction through solids and by radiation and convection through air spaces. 
In addition, the program has built-in features for ease of application to 
building constructions where various combinations of solid-to-solid and solid- 
to-air interfaces are encountered. For each solid layer, the present program 
has provisions for phase changes, heat generation and absorption, and thermal 
property variation with temperature commonly found in building materials. 
Through the air spaces the modes of heat transfer include radiation and con- 


vection with temperature-dependent air properties. 


A number of special analog and numerical programs for fire endurance 
evaluation have been in existence for some time fisay.. A more flexible and 
general finite difference program was developed by Krokosky as recently as 
1970 [5], which also contains a review of the existing thermal analysis 


programs for fire endurance evaluation. 


The present numerical program was also developed to incorporate into 
fire endurance analyses the following features which are desirable and yet 


not readily available in existing programs: 


1. Options to handle heat exposure on one or both sides of a structure. 

2. Allowance for air infiltration, and heat generation and absorption 
in air spaces. 

3. For ease of application to building structures, a single input card, 
using digits 1 and 0, is sufficient to instruct the computer the 
specified number of solid layers and air spaces arranged in a given 
problem. 

4, Temperature-dependent thermal properties and known heat effects due 
to chemical reactions of some common building materials are stored 
in a subroutine and called by an input card, punched as 2331, 
indicate coded materials stored in the subroutine. 

5. The duration as well as the amount of heat represented by a phase 


change or chemical reaction can be varied in any material. 


t Numbers in brackets refer to the literature references listed at the end 


of this paper. 


2. GOVERNING EQUATION AND BOUNDARY CONDITIONS 


The governing equation for one-dimensional transient heat flow is the 
well known heat diffusion equation. Including a term for internal heat 


generation, this can be expressed as: 


Se eee a 

where 

T = absolute temperature in solid 

X = coordinate in direction of heat flow 

k = heat conduction coefficient 

g = time rate of heat generation per unit volume in solid 

oe = density of solid 

c = specific heat capacity of solid 

6 = time. 


Boundary Conditions [6,7]: 


The following types of boundary conditions need to be evaluated for the 


general problem: 


(ayresolid to air 

(b) energy balance in air space 
(c) air to solid 

(ad) solid to surrounding 

(e) solid to symmetry plane 

(2) 4sollid-‘torsolid 


(g) furnace gases to first surface layer. 


AX, AX, 
< Lar 
zi Ta 79 
Air 
Solid 1 Space Solid 2 


Figure 1. Solid-Air-Solid Interface 


(a) ~ Solid to: Arei(see (Eagm ih) sail 


—- = ke ne =i'gle, G = =) S Biya Ds a (2) 
(b) Energy Balance in Air Space (see fig. 1) 

NaCl oe n(?,-7))°7* = n(n ath) 7 om co (rn aie (3) 
(c) Air to Solid (see fig. 1) 


In the above equations the subscripts 1 and 2 indicate solid 1 and solid 2, 


respectively, as shown in figure 1, and where 


ip 
Slate ees ar 
St oi 
and 
ey E, are emissivity of surface 1 and of surface 2 
T,, Ty? surface temperatures as shown in figure 1 
AX), AX, : incremental length in x direction in solid 1 and solid 2 
oO: Stefan Boltzmann constant 
ne: temperature of air 
a, eee heat conduction coefficient in solid 1, and 2 
Gar ede time rate of heat generation or absorption per unit volume in 
solid 1 and 2 
hi air gap spacing 
a: specific heat capacity of air 
Po: density of air 


m: rate of mass transfer per unit area due to pressure difference 


go: rate of heat generation per unit volume in air space due to 
combustibles 

h: an empirical convection heat transfer constant, 
h = .27 for vertical surface, and h = .38 for horizontal 


surfaces [8]. 


A 
Bedsaeton 
Solid T A 
——— 
Convection 
Figure 2. Solid-to-Surrounding Interface 
(d) Solid to Ambient (see fig. 2) 
ocAx De Oa AOS Cs se x 5/4 Ax 
mon. AG rs oe(T T,) hie fe) +95 (5) 
where the additional variables are: 
Bos ambient temperature 
E: emissivity of surface 
aus surface temperature of solid 


one” plane 


= AS = Ax 
2. 
Solid 
Solid 
nee 


Figure 3. Symmetry Plane 


(e) Solid to Symmetry Plane (see fig. 3) 


By symmetry, temperature on both solid surfaces that face each other are 
equal, however T > To. So there will be heat transfer from solid to air. At 


the interface we have: 


DCA ATs eee C ee " 5/4 Ax 
Oa AGT ait k ae h(T TS) ce Ones (6a) 
At the air space we have, 
AT 
5/4 
© eho. = Sener 4 (6b) 
S paa A@ s 


where 
Es surface temperature of solids 
TS: temperature of air in symmetry plane 
he: distance from solid to symmetry plane 
TW 22 
Solid 1 Solid 2 


Figure 4. Solid-to-Solid Interface 


(£) Solid to Solid (see fig. 4) 


Consider interface between Solid 1 and Solid 2 as shown in figure 4. 
Let T, and T, be the interface temperature in each solid respectively. A 


heat balance for the interface node can be expressed as u = T , and 
1 


Ax AT, AX, AT, Ax, Ax AT AT 
pe 2 ee onc a 0) =e ee a Se oe eee k, — 


2 NG 149 et 


(7) 


where subscript indicates conditions in solid layer 1 or 2. 


Radiation 
Th qT) Solid 
Convection 


Figure 5. Furnace-to-First-Solid Interface 
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(g)i Furnace Gases to First Solid Surface (see fig. 5) 


Consider heat transfer from furnace gases to first solid surface. The 
main modes of heat transfer will be radiation and convection. A good approxi- 


mation of ASTM E 119 fire curve is given by the following three formulas: 


Toe= W692 == + 68 Oj* 1° <*50min (error +. 2%) 
ie 1667 +. 1.260 - 0.0236 (120-6). 50 min, < 68° <\ 115° min 

(no appreciable error) 
Tae l667 tel «260 115 min < 6 < 480 min 


(linear exact) 


Heat balance from furnace to first solid surface, 


pcAx : = k AT 
7 ING) Ax 


. + 5/4 AX 
a = Lis 
+ oe (T,, ao a h(T, oe arc ae (8) 


where Tr is absolute furnace temperature. 


3. PROBLEM FORMULTION AND FINITE DIFFERENCE EQUATIONS 


Consider general one-dimensional heat transfer problem containing N 
solid layers and m air spaces in any order. To facilitate formulation and 


discussion we introduce i, j, and k, three running indices. 


Let: 
ae: Temperature of jth node in ith solid 
Te: 1...N (number of solid layers) 
ae 1...n (number of nodes in solid layer 
ks 1...m (number of air spaces) 


In figure 6, a general multi-layer configuration is shown, where numbers in 
circles indicate the applicable equation at the given node as discussed in 


the previous section. 


Radiation 


Convection 


Figure 6. A Composite Wall with N Layers, m Air Spaces 
(Numbers in circles indicate applicable equation 
as discussed in section 2.) 
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To solve our problem numerically with the nonlinear heat diffusion equation 
and associated complex system of boundary conditions we shall use forward 
time differencing and central space differencing scheme. Figure 7 shows the 


finite differencing network inside a solid layer. 


CG 
~— 


nl 


i 


| 

I 

| Ax 

2 nth node 
| b 

| 


| 
| 
\ 
Pine 
1 
j=l...n 


ith solid layer 


Figure 7. Sketch of Finite Differencing Network in a Solid Layer 


Furthermore, let 


Tok? Temperature of kth air spacing 


6: Time 
A868: Time increment 
Ts ne Temperature at (6 + A686) of jth node in ith solid 
, 
Ras Thickness of ith solid 
dk? Thickness of kth air spacing 
Oe 
Ax, = + : Difference spacing for ith solid 
gs os Rate of heat generation per unit volume 
v 
k, a Heat conduction coefficient 
v 
Os Heat diffusion coefficient 


Sak? Rate of heat generation per unit volume in air space 


fe} 
my)? Rate of mass transfer per unit area 


Applying our finite differencing scheme we have the following general 
finite difference expressions corresponding to the previously discussed 


governing equation and boundary conditions [6,7 and 9]. 


Governing Equation 


2 
ak &: 2 g, Ax, 
Te Si oT. fA ral 8%, oo a cae | aa 
ah M, ( Li(ai=L) ie) sh ( x) Ti5 + kom 3} 
i Dok 
where et ) 2 
ee a) 
ioe OMA 
dl: 
From equation (9) we require M; > 2 for stability. 
Solid to Air 
rh = T $28 nL - T tae cL; : T ; 10 
in Spi Ma (Td a1 Sal i 1 Ty (Toya) eu 
: A 
- H (T, T,)°/* 4 stun Ti 
Sey elec) ak : ‘ 
i ota 
where 
i é 2AGce,, 
ts ponremenn 
ig id: 
rr: Specific heat of ith layer 
2A6h 
a (b,c, Ax, 
Agar Ga 
Air Space Heat Balance 
5/4 5/4 + G 
' — — _ 
wai) crak °° Ti,n ak! a Serer atmo a or a eed eS 
CE) 
where 
hA8 
P =) 
‘ "axopa a 
Po: air density function of temperature 


GJ: “specific heat of air 


pa 
hA@ 
(@) = 
> ‘ak°pa’a 
fe} 
mA 9 
S = 
Qg 
k aka 
A8 
G = Jak 
ak Chap 
paa 
Air to Solid 
T! = T, ae ea | CidPR pee (el Yas ete (a ~ 
it+1,1 itl ,1 ne in ae bases M, sista al 
itl 
< 5/4 
a ee ee ne 
Rey Misi, Head as defined before. 
Solid to Ambient 
T ST eT = pe. SRG ene yarn 
on N,n M, N,n-1L N,n N voNiy,n NON 
2 
A 
ee Xyy) 
+ 
nN 
where SNelate 
N 
Ry = —— 
AX Pon 


emissivity of solid (last node to ambient) 


N° 
Hy Nn and My as defined before 
TO? ambient temperature 
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(12) 


(13) 


Solid to Symmetry Plane 


AX,.) 
2 5/4 n,n (o%y 
an = edit a (Pet ape yy Se aT = 7) oo 
N,n N,n My N,n-l N,n N° N,n s npnn 
a mes 4 
PNoNen ots} (14a) 
Bs 5/4 4 4 
' = + — - 
T ub Ply aS) + Ry (Ty on T, ) (14b) 
where 
Ta temperature at symmetry plane 
A8 
a Te 2 0 
Ss-pa’a 
Xa: distance between solid surface and symmetry plane. 
Solid to Solid 
' = ' 
Gee ign 
1 — - 
ieee 7 fae eet eri on om eee eg ey Daeg Ee ene Edie) 
oe aia gae eS a are Aa sy Bey) (12) 
where 
i 1! pc, Ax, 
al ay 2 
Ax. A8® 
Bityes == 
z 2 
kj ne 
Dein ASK. 
i 
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Gia ais rate of heat generation per unit volume in ith solid. 


Furnace Gases to First Solid Surface 


E i Cities 
] = a y 
Dy ag Re ee ra ec oT 
allie AL 
~ 220 a ee | 16 
M, 1,1 12 (16) 
where 
248cc. 
asl ie x Cc 
eee y ds 
ee emissivity of first surface (furnace gases to lst solid). 
Th furnace temperature, and 
Hy: as defined before 


4. DISCUSSION OF NUMERICAL PROGRAM 


A flow chart of the main program is presented in figure 16. A listing of 
the complete numerical program is also presented at the end of this report. 
The logical sequence of the main program can be grasped readily by first 


considering the flow chart. 


The following is a list of input parameters as read in on data cards 


(examples of which are shown at section 5.3). 


NN: number of solid layers 
N: number of nodes in each layer 
M: number of air spaces 


ID: sequence of numbers 1 or 0, specifying the sequence of solid layers 
and air spaces, e.g., 101011101 means solid-air-solid-air-solid-solid- 
solid-air-solid. 
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IDD: sequence of numbers specifying the material of each solid layer, 
e.g., 122321 means the six solid layers of the problem are of 
Materials 1, 2, 2, 3, 2, 1 in that order. Various temperature 
dependent material properties are stored in Subroutine PROP. 

AL(i): thickness of ith solid layer in feet 

AI(i): thickness of air gap spacing in feet 

GA(k): rate of heat generation per unit volume due to combustibles in kth 
. air space, Btu/ft?-hr. 

DTHETA: Time increment in fraction of hour. 

AMO: Volume rate of airflow through walls per unit area in ft°*/ft?-hr 
(m/p) . 

GM1, GM2, GM3, etc: Heat absorption or release due to phase change in 
Matcecialy 155 2Qvort 3.etc,.,. Btu/ft*hr. 


Subroutine PROP stores temperature dependent thermal properties and phase 
change reactions of some common building materials. This subroutine can be 
expanded as desired when new materials are encountered. When called from the 
main program this subroutine supplies the thermal properties for the ith layer 


of solid currently being calculated. 

The function subroutines are self-explanatory. The comment statement at 
the beginning of each function subroutine properly identifies it with the 
corresponding equation and boundary conditions in section 3. 

5. SAMPLE CASES CALCULATED 


5.1. Description of Panels 


Case 1: Four solid layers (gypsum, wood, wood, gypsum), three air gaps, heat 


from one side as shown, 


where 
G: Gypsum board 
W: Plywood 
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Case 2: Four solid layers (gypsum, wood, wood, gypsum), three air gaps, 


heating from both sides as shown, 


where 


G: Gypsum board 
W: Plywood 


Case 3: Bulkhead 2B”, four solid layers (marinite, mineral wool, aluminum, 


marinite), two air gaps, 


heat from one side as shown, 


M : Marinite 
MW: Mineral Wool 


AL: Aluminum 


Case 4: Bulkhead 3, three solid layers (marinite, aluminum, marinite), two 


air gaps, heat from one side as shown, 


where 
M : Marinite 


AL: Aluminum 


Assembly descriptions and details of the actual ASTM E 119 fire endurance 


tests of this and following bulkheads are presented in reference [11]. 
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marinite), two air gaps, heat from one side as shown, 


| - Air 


6" 


M :. Marinite 


AL: Aluminum 
MW: Mineral Wool 


Case 5: Bulkhead 3C, four solid layers (marinate, aluminum, mineral wool, 


Case 6: Bulkhead 3F, four solid layers (marinite, mineral wool, aluminum, 


Marinite), two air gaps, heat from one side as shown, 


M 


Case 7: Bulkhead 4, three solid layers (perforated steel, mineral wool, 


aluminum), one air gap, heat from one side as shown, 


Air 


a 


where 


PS: Perforated steel sheet (#16 USCG, 3/16 inch diameter holes, 
3/8 inch diameter holes, 3/8 inch center, 23% open area) 
MW: Mineral wool 


AL: Aluminum 


Case 8: Deck 5Bl, three solid layers (perforated steel, mineral wool, 


aluminum), heat from bottom as shown, 


Air 9" 


PS 


AB's) 


5.2. Thermal Properties 


The following thermal properties are taken from Ref. [1] and Ref. [10]. 


NBS property measurements are programmed into subroutine PROPT, 


Gypsum Board 


p = 60 lb/ft” 

c = .26 Btu/lb °F 
= .l25, Btu/hr, ft, 38 0) < T s 20000r 
= .075 Btu/hr ft °F 200 °F < T < 400 °F 
= (.05 + Te7G00) Btu/hr ft °F 400 °F < T < 2000 °F 

eee ied) 


Heat of moisture desorption and calcination: 20,740 Btu/ft 


Brick 
o = “Lio Tbr 
ec = 216 Btu/lb- SF 
eee 9 
k = 1.0 Btu/hr ft °F 0 <P < | 2004eF 
= (.46 42771:0;:000): Btu/hr Eo SE 200 °F < f < 2000 
Heat of moisture desorption 5000 Btu/ft? 
Plywood 
6 =" 31 ab/rtte 
© =) soos Btu/ ib ck 
—e = .9 
k = 0652 Btu/he hts or 
Heat of moisture desorption 1500 Btu/ft? 
Mineral Wool 
p = 6.04 lb/ft? 
GC = .2°>Beu/1lb. or 
e = .9 
k = 02> Btu/hre see 0S) Tes, 510) 
= 9.02167 Beu/hrest wor SUF < <7 Sack 
= .0229 + 0.517 (T-100)/10,000 (SESE s T <<) 2000R ee 
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aa 


Aluminum 


se338 f£b2/hr 


Q 
Il 


pie Sta 69, Wib/ Et 3 

CH= a 2UShBEU/ Lb: OF 

cae) 10/9 

em. Loy. + 1.63 T/10, 000 0 < "Tax 200), °F 
Se a LS (T= 32) 200 << 2000; OF 

Marinite 

—we36. Tb/ ft? 

Coen .25eBtu/ib) °F Ov <“200  °°F 
zor: 3(l->200)/10;,000 Z000C°F -< =<" 20007°F 

eee =H 9 

koe—oeeOG3 Btu/hre ft °F so <= - 200 TF 
—meniea7io to. 1 ob (T-1L00)/10,000) /12 200 °F <"T-< "2000 °F 


Saas “Input: Data 


The following are the complete input data for the sample cases calculated. 


Case 1: Four solid layers, 3 air gaps, furnace on one side (gypsum-air-wood- 


air-wood-air-gypsum) 


4531 68. -0 513:8.0)...0 1500.0 310 
1010101 

33: 

= (54 .0417 -0417 7052 

a2) O. evs) OX. BPA) Oi. 


Case 2: Four solid layers, 3 air gaps, furnace on both sides (gypsum-air- 


wood-air-wood-air-gypsum) 


4532 68.0 .0 5380.0 1500.0 185 (0) 
1010101 

1g )1 

OD)Z -0417 -0417 -052 

AS) 0. x22 0. 0. 


17, 


Case 3: Bulkhead 2B, four solid layers, 2 air gaps, furnace on one side 


(marinite-air-mineral wool-aluminum-air-marinite) 


4321 68.0 -0 -0 5 (0) IE) 

LOLLOL 

6456 

063 .083 0.052 083 

25 .00 . 083 .00 

Case 4: Bulkhead 3, three solid layers, 2 air gaps, furnace on one side 


(marinite-air-aluminum-air-marinite) 


3321 68.0 .0 HO .0 10 
10101 

656 

.063 .052 063 

083 0 5 .0 


Case 5: Bulkhead 3C, four solid layers, two air gaps furnace on one side 


(marinite-air-aluminum-mineral wool-air-marinite) 


4331 68 -00 -00 -00 1310 
101101 

6546 

063 -05 -083 042 

083 -00 at) -00 


Case 6: Bulkhead 3F, four solid layers, two air gaps, furnace on one side 


(marinite-air-mineral wool-aluminum-air-marinite) 


4321 68.0 -00 -00 .00 130 
101100 

6456 

-063 - 083 «052 063 

083 -00 i) - 00 
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6. DISCUSSION OF RESULTS AND CONCLUSIONS 


Comparison plots of calculated and measured temperatures from standard 


fire endurance tests represented by the above sample cases are presented in 


figures 8-15. Tests were conducted at the NBS large scale fire test center, 
Washington, D.C. Complete specification of the bulkheads and decks and 

descriptions of the tests are presented in a separate NBS report [11]. The 
good agreement between calculations and actual fire endurance test measure- 
ments as shown in the comparison plots reflect the success of the mathemati- 


cal modeling. 


In the comparison plots, temperature-time histories were recorded for 
each individual panel in a given wall construction for comparison with 
computer calculations. The plots can be examined for details of the heat 
transfer processes. For the purpose of fire endurance rating generally only 
the temperature history of the unexposed surface is needed. Table 1 shows a 
comparison of the fire endurance times by experiment and calculation. The 
criterion used in this table is the 250 °F average temperature rise condition 
at the unexposed surface. The agreement of experimental results and theore- 
tical predictions of fire endurance times is good. Based on those calcula- 
tions one would recommend using the present computer program for the purposes 
of fire endurance evaluation and design studies. Such applications would 
result in considerable savings in both money and time. This program could 


also be used to evaluate the performance of new materials. 


Finally a word of caution in the use of a thermal analyzer program 
with air gap heat balances such as the present program. Due to the low heat 
Capacity of air, one must exercise caution in selecting the incremental Aé. 
The iteration time, AO, used in the previous sample calculations are 1/720 or 
1/1440 hour for cases which do not involve metal layers and 1/14400 hour for 


cases involving metals. 
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Table l. 


Comparison of Fire Endurance Times by Experiment and Calculations 


Description 


4 layers, Gypsum, Wood, Wood, 
Gypsum, 3 Air Gaps 


4 Layers, Gypsum, Wood, Wood, 
Gypsum, 3 Air Gaps, fire on both 


sides 


4 Layers, Marinite, Mineral Wool 


Aluminum, Marinite, 2 Air Gaps 


3 Layers, Marinite, Aluminum, 


Marinite, 2 Air Gaps 


4 Layers, Marinite, Aluminum, 
Mineral Wool, marinite, 2 Air 


Gaps 


4 Layers, Marinite, Mineral Wool, 


Aluminum, Marinite, 2 Air Gaps 


3 Layers, Perforated Steel, 


Mineral Wool, Aluminum 1 Air Gap 
3 Layers, Perforated Steel, 


Mineral Wool, Aluminum, 1 Air 


Gap 


20 


Endurance Time 


Passed 1 hr 


Passed 


Passed 


Passed 


Passed 


Failed 


38 min 


Failed 


30 min 


calculation 


Passed 1 hr 


Failed 


37 min 


Passed 


Passed 


Failed 


36 min 


Failed 


34 min 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


baal 


[8] 


[9] 


[10] 


fell) 
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Figure 12. Thermal Response of Bulkhead 3c. 
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Figure 13. Thermal Response of Bulkhead 3F. 
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